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We investigate the effect of electron- Ai 9 phonon coupling on the gapless electronic band dispersion 
of the pristine graphene. The electron-phonon interaction is introduced through a Kekule-type 
distortion giving rise to inter-valley scattering between K and K' points in graphene. We develop a 
Frolich type Hamiltonian within the continuum model in the long wave length limit. By presenting a 
fully theoretical analysis, we show that the interaction of charge carriers with the highest frequency 
zone-boundary phonon mode of Ai 9 -symmetry induces a mini band gap at the corners of the two- 
dimensional Brillouin zone of the graphene. Since electron-electron interactions favor this type of 
lattice distortion, it is expected to be enhanced, and thus its quantitative implications might be 
measurable in graphene. 

PACS numbers: 71.38.-k,63.22.Rc,72.80.Vp,78.30.Na 



I. INTRODUCTION 

Since the discovery of graphenei^ and its unconven- 
tional physical properties, the investigation of electronic 
properties of both graphene and graphene based nanos- 
tructures have become one of the active areas in con- 
densed matter physics experimentally as well as theo- 
reticallyin past few years. Theoretically, in the low- 
energy limit charge carries of graphene have linear disper- 
sion relation around so-called Dirac points 3 - having Fermi 
velocity- vp ~ 10 6 m/s, and Dirac- Weyl equation can be 
safely used within the framework of continuum descrip- 
tion of the electronic band structure of the graphene^. 

It is also well-known that both in-plane and out-of- 
plane phonon modes play an important role in charge 
carriers dynamics of the graphene^—. On the one hand 
the Fermi velocity is reduced by interaction of charge 
carriers with doubly degenerate in-plane Fi2 g phono n 27 ' 28 . 
On the other hand, though the electron-highest frequency 
zone-boundary phonon interaction, i.e., Kekule-type dis- 
tortion of the graphene lattice is one of the possible mech- 
anisms among the gap generations, except that the work 
of Reffl6[, there are no theoretical works about its influ- 
ence on the graphene band dispersion. This first theo- 
retical prediction of dynamical mini band gap formation 
in graphene due to the highest frequency phonon mode 
with Ai s -symmetry is reported by Samsonidze et alJ£. 
They showed that, based on a simple tight-binding model 
at room temperature, such as electron-phonon coupling 
mechanism induces a mini gap around 10 meV, and it is 
also responsible for the Kohn anomalie s 32 ' 33 in graphene. 
The Kekule structure consists of a network of hexagons 
with the alternating short and long bonds like in the 
classical benzene molecule. This pattern was studied for 
ID simple model, finite size carbon nanotube a 34 ' 35 . In- 
vestigation of the gap formation, in particular, its con- 
trol, in both graphene and graphene based nanostruc- 
tures is itself one of the hot topics of the current research 
in graphene, and such a gap generation can be created 
by strain 3 ^— or by substrate induced effect o 10 ' 39 ' 40 . 

In this paper, to investigate the effect of interaction of 
graphene charge carriers with the highest frequency opti- 



cal phonon mode of Ai ff symmetry near the zone bound- 
ary K (K'), we performed an analytical study based on 
Lee-Low and Pines (LLP) theory 4 ^.. The carrier-phonon 
interaction is described through a Kekule-type distortion 
giving rise to inter-valley scattering between K and K' 
points in graphene 4 ^— . Based on this interaction, we 
first construct a Frolich type Dirac- Weyl Hamiltonian 
which is nondiagonal in phonon creation and annihila- 
tion operators. Secondly, we present a simple analytical 
model to diagonalize it by just introducing two succes- 
sive unitary transformations. Finally, we show that the 
interaction of charge carriers with highest zone boundary 
phonon mode opens a mini band gap at the corners of 
the Brillouin zone. 



II. THEORY 

In the long- wave length regime, the Hamiltonian of 
the graphene electron (hole) interacting with Ai s -phonon 
mode can be written as 
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where = vpot-p is the unperturbed part, whose spec- 
trum describes cone like behavior around the Dirac points 
with eigenvalues — Xvpk. A is the chirality index, and 
takes — 1(+1) values corresponding to valence (conduc- 
tion) bands in pristine graphene. These two bands touch 
each other at the corners of the Brillouin zone, i.e., at the 
well-known K and K' points whose coordinates are given 
by K = (2n/a) (l/3, and, K'= (2tt/o) (2/3,0), re- 
spectively. We have labeled these points in Eq.(JT|) by 
the valley index ji. Here, a are the four component 
Dirac matrices, and a is the equilibrium bond length, 
i.e., 1.42 A. Thus, the corresponding eigenfunctions of 
the unperturbed part can easily be constructed in 
terms of four component pseudospinors 
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where I? is the total area of the system. In Eq. (fTJ) , the 
last term represents the electron- phonon couplings 45 , 
and is given by 
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where /?k = — din Jo/<ilna, 7 = (3a/2)Jo, ui — 
exp(2-7ri/3), Jo is the resonance integral between near- 
est neighbor carbon atoms which is of order of 2.77 eV, 
dy is the 2x2 Pauli matrix. In Eq.Q, the amplitude of 
distortions at K and K' points are defined by 
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respectively. N is the number of unit cells, Mc is the 
mass of a carbon atom. In Eq.Q, 6k. q (&K',g) an d 
b K q (^k' q) are the phonon creation and annihilation 
operators at points K (K 7 ) with phonon wave vector q 
and frequency uik{q)- The corresponding highest zone- 
boundary phonon energy is 7Lo>k:(0) = 161.2 meV. There- 
fore, the electron-phonon interaction Hamiltonian given 
by Eq.(|3]) can be conveniently rewritten in the following 
form: 
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Here, ao = 
(dJ /da)/J is 



together with M = 3a <7o<A) 
(^/2M cWK (0)) 1/2 , and q = 
predictedil^, around 2.0 A -1 and 2.5 A' 1 . 

To diagonalize the phonon subsystem of Eq.([T]) 
through Eq.Q we employ a unitary transformation 



scheme within the LLP theory. This includes two succes- 
sive transformations each of which eliminates the electron 
coordinates from Eq.(fT]), and shifts phonon coordinates 
by an amount of the interaction strength, respectively. 
To do this we follow the method developed for the inves- 
tigation of the interaction of electron (hole) with doubly 
degenerate optical phonon modes of E2 9 symmetry near 
the zone center—, wherein an ansatz was made so as to 
take into account chiral nature of the pristine graphene 
due to its gapless electronic band structure. However, 
besides the chiral nature of the problem, one must also 
considered that the zone boundary phonon gives rise to 
inter-valley scattering between K and K'. Therefore, to 
be compatible with these properties of the problem, we 
make an ansatz for the ground-state of the whole system. 
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such that H \ <&) = E± | <&). Here, | 0) p h stands for the 

phonon vacuum, and a± A | fi'Xk) corresponds to elec- 
tronic state vector defined through the appropriate frac- 
tional amplitudes, a±*, due to the fact that total wave 
function of the system must be the linear combination of 
I fi' + k) and | p! — k) , respectively. 

On the one hand, the first unitary transformation 
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eliminates electron coordinates from Eq.([T]), since the 
transformed operators are given by the relations, 6 Mi q = 
b^g exp [-iq • r] and p = p-J2J2 ^Qb^qb^q. There- 

q njtu 

fore, the transformed Hamiltonian takes the form, 
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On the other hand, second unitary transformation 
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is the well-known displaced oscillator transformation 
which shifts phonon coordinates by an amount of the 
interaction amplitude, Mo =Mo/7Lu;k(0). It just shifts 
the phonon coordinates, since it generates the coher- 
ent states for the phonon subsystem such that opti- 
cal phonon operators transform according to the rule 
bfi.q = bn >q + M (fi'X k I ML I v'Xk). As a result, un- 
der the transformation U2, Eq.® can then be written as 
H = H° + Hi, where H° and Hi are given by 
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respectively. Therefore, one applies the phonon vacuum 
to the sum of Eq. (fT0)) and Eq. (fTTj) . only the contribution 
comes from the diagonalized part , i.e. from H. . By us- 
ing the ansatz given by Eq.© , one first applies Eq. (fTU|) 

fj, A 



struct the eigenvalue equation % | 3>)^ A = E± \ <^)i xuX . 
Finally, by taking inner products to compare the related 
coefficients of the states | //A'fc) we arrive four simultane- 
ous equations for , a!j5~ 
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to the term a± | //A'fc}, and then sums over A' to con- be rewritten in the following matrix equation: 
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where we have defined Oaa' as 



@AA' = - [(« + s') cos(e -<p)- (s - 
together with the matrix elements 
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where we have defined that ao — \Mq\ /4 JoHujk(0) which 
takes values 0.305 or 0.477 depending on whether q Q is 
2.0 A -1 or 2.5 A -1 , respectivel y 47 ' 48 . We must introduce 
an upper cut-off frequency q c = q c a, while taking the in- 
tegrals in Eq. (|13[) . since they diverge at upper limit of the 
integrations. By solving the determinant of the matrix 
in Eq.(fT2"| the eigenvalues E± can be solved analytically 
in the following form: 



which are equal to unity for C 7^ C an d A ^ A', other- 
wise zero. This shows that only inter- valley scattering 
having different chiralities are allowed due to the con- 
servation of the chiral symmetry . In Eq. (fl"4|) . 9(<j>) are 
the azimuthal angle of the momentum k{q). After con- 
verting the sums in Eq. (|13[) into integrals over q, i.e., 
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which is modified electronic band dispersion of the pris- 
tine graphene by the gap due to the Kekule-type distor- 
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tion of the lattice. As is seen from Eq. pTl) . Kekule- 
type distortion preserves the chirality of the sublattice, 
i.e., the valley degeneracy is not lifted. Since the calcu- 
lations we have done are restricted by the energy scale 
near to the phonon resonance, i.e., wk, we can choose 
q c = uik/vf = 0.027 A" 1 such that q c = q c a = 0.039. 
This suggest that the magnitude of the half-band gap 
is of order of 2.12 meV (q = 2.0 I" 1 ) or 3.34 meV 
(go = 2.5 A ). so that the induced gap, of 4.24 meV 
or 6.68 meV. This is smaller than that those previously 
found in the literature^, where a mini band gap occurs 
10 meV in their room temperature calculations. 

III. CONCLUSION 

In conclusion, we present, to the best our knowledge, 
the first theoretical justification of a mini band gap for- 
mation in pristine graphene at absolute zero tempera- 
ture, due to the interaction of electron (hole) with high- 
est frequency optical phonon mode with K\ g symmetry 
near the zone boundary yielding Kekule type distortion. 
We have shown that such an interaction opens a gap 
without breaking the chiral symmetry of the lattice. Al- 
though, we have performed our calculations at absolute 



zero, it is compatible with those found at room tempera- 
tures for both graphene and zig-zag single walled carbon 
nanotubes^. In the graphene literature, there are only 
two theoretical work o 46 i 49 and no experimental work de- 
voted to the investigation of combined effects of electron- 
electron interaction and electron-Ai s phonon interaction. 
On the one hand, since the electron repulsion breaks the 
sublattice symmetry whereas the Kekule patterned lat- 
tice distortion restores it, it is predicted in Ref. H^that , 
taking such competing interactions, superconducting or- 
der maybe have to be considered. On the other hand, 
Giuliani et al4^ predicted that the electron repulsion en- 
hances dramatically, the gaps due to the Kekule distor- 
tion. We believe that the results of this paper may be 
useful for understanding the role of electron-electron in- 
teractions in graphene. 
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